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Abstract 

' We show that the manifold of density matrices can be derived from 

Cn . CP^ by the action of SU{N). We give some prehminary observations 

on the structure of this manifold. 
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1 Introduction 



The manifold of pure quantum states is well known and well described, but 
the generalized quantum states, the set of density matrices (hermitian positive 
matrices of unit trace) remain not very well understood. 

In this paper we set out to find the manifold of density matrices. We show 
that the manifold of density matrices can be derived from the CP^ -1 by the 
action oiSU(N). 

2 Manifold of density matrices 

The technique of purification allows us to describe NxN density matrices in 
terms of dimensional pure states. Any NxN density matrix has a canonical 
decomposition: 

N 

p = (1) 

i 

where pi and are the eigenvalues and eigenvectors, respectively. 
Using the purification procedure, we can always find a pure state in a N'^ 
dimensional Hilbert space, given by its Schmidt decomposition: 

N 
i 

where the partial trace over the ancillary system yields the original density 
matrix: 

N 
i 

Note that we have labelled the states in the original space of the density 
matrix with the superscript A and the ancillary system with the superscript B. 
Also note that forms an orthornormal basis. 

We can consider this as a many-to-one mapping from the set of pure states 
in N'^ dimensions to density matrices in N dimensions. The following lemma 
allows us to determine just how many pure states would correspond to the same 
density matrix. 

Lemma 1: Given 2 pure states in a N'^ dimensional Hilbert space \(f>) and 
lip), with the same partial trace: 

Trsmm^TrBrnm (4) 
then there exists a unitary transformation v e SU{N) such that /"^^w"^ 1-0) = 

10)- 
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To prove Lemma 1, let us write the 2 states in their Schmidt decomposition: 

l^> = E*k,^>l9f) (5) 

k 

\4>) = jZpk\piM) (6) 

k 

Their partial traces are equal: 

Y,<il\<it){qt\ = T,pl\pt){pt\ (7) 

fc k 

iS — i^, ql — p1 and if no values ql are degenerate then \q^) = \p^) upto 
a phase. 

In the non-degenerate case, the only freedom is in speciiying the basis {\q^)} 
and )}, which is given by a unitary transformation in the B system. 

If a degeneracy exists between for particular values of z, then the states 
\q^) and \p^) can differ by an unitary transform. We will show that this unitary 
transformation can be equivalently performed in the B subspace. Consider the 
subspace of Itjj) for which the Schmidt coefBcients are degenerate: 

(8) 

z 

A unitary transformation R within the space spanned by \q^) gives: 



Y,CR\qM) = E^E^-'l9.^)laf) = T.CT.\l^)R.'Mf') = 

z z z' z z' z 

(9) 

Therefore we see that i? (g) /|^) = / ® 

And finally to complete the proof, any complex phase on the Schmidt coef- 
ficients can also be introduced by a unitary transformation on the B subspace. 

With lemma 1, we see that SU{N) orbits within the manifold of N"^ dimen- 
sional pure states CP^ ~^ correspond to unique N x N density matrices. This 
is a bijection between the set of N x N density matrices and these orbits on the 
manifold CP^^-^/SU{N). 

3 Geometry of the manifold of density matrices 

Let us consider the action of SU (N) on the manifold CP^^~^ using the follow- 
ing: 

I^<^U^Y.qk\qtM) (10) 
k 
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where S SU{N) acts on a state in CP^ ~^ written in its Schmidt 
decomposition. 

This action is not free, therefore we need to classify the stabihzers of the 

states. For a state; with Schmidt mimbcr //,, the orthornormal states \qj^) would 
span a /U dimensional subspace, so the stabilizer of this state is the subgroup 
U{N — fj). This gives rise to different dimensional strata, identified by the 
Schmidt mimber (or equivalently the density matrix rank) /i. 

The dimension of CP^'-^ is 2N'^ - 2 and the dimension of SU{N) is iV^ - 1. 
This gives us the expected dimension of N"^ — 1 for the density matrices. For 
states with maximal Schmidt number ^ = N , the stabilizer is the trivial group 
1, the dimension of the stratum is A^^ — 1. As the Schmidt number decreases, the 
stabilizer group increases, and the dimension of the stratum decreases. For rank 
matrices, the stratum is ji{2N — /i) — 1 dimensional. The lowest dimensional 
stratum is 2N — 2, which is the set of pure states, and which is also the manifold 

Each stratum is also a convex covering for the higher dimensional stratum, 
as it is clear that a rank /z density matrix can always be decomposed in terms 
of rank /x — 1 density matrices. The lowest stratum, of the set of pure states, is 
the overall convex covering for all density matrices. 

We now see the set of density matrices as a layering/strata of surfaces of 
dimensions given by ii{2N — fjb) — l, with 1> li> N , with the lowest dimensional 
stratum, a 2A'' — 2 dimensional surface forming the set of pure states. 

4 Geometry of N=2 density matrices 

As a simple example, let us consider the case where N = 2, the manifold of 
density matrices is given by CP^ / SU(2). There are only 2 possible density 
matrix ranks, giving rise to only 2 strata. The stratum of pure states forms the 
convex hull and is 2 dimensional. The stratum of rank 2 matrices, the mixed 
states, is 3 dimensional. 

The manifold is symmetric under 5/7(2), which coincides with 5*0(3) x 
Z2- Given the symmetry and structure, the only possible shape is a sphere. 
This agrees with the well known Bloch sphere representation for qubit density 
matrices. 

5 Conclusions and comments 

In this paper, we have obtained an interesting relationship between the mani- 
fold of density matrices and the well known manifolds of pure states CP^ . This 
allows us to discuss a mathematical manifold for density matrices which has not 
been available to us before. With this, we hope to develop a greater under- 
standing of quantum density matrices and their operators (dynamical maps). 
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